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Abstract  -  In  the  resonance  theory,  the  radial  dilatation  is 
emphasized,  and  the  blood  pressure  wave  is  transmitted  in  the 
form  of  “moving  windkessel”.  Based  on  this  conjecture,  we 
developed  a  semi-empirical  procedure  to  describe  the  pressure 
distribution  in  a  complex  simulated  model  composed  of  a  main 
tube  and  attached  organs.  Now  we  try  to  apply  this  fitting 
method  to  hemodynamics  in  vivo,  and  we  tested  our  equations 
by  two  sets  of  experiments  on  the  rats:  ligating  a  renal  artery  and 
bending  the  aorta  transversely.  Abdominal  aortic  blood  pressure  of 
rat  was  measured  through  tubes  inserted  fom  the  caudate  artery 
and  compared  with  curve -fitting  deduced  from  semi-empirical 
resonance  equations.  The  good  fitting  result  illustrates  that  in  spite 
of  various  complex  structures  of  the  arterial  system,  we  can  still 
provide  good  description  for  the  blood  pressure  distribution  by 
dividing  the  arterial  system  into  sub-units  and  describing  with  few 
elastic  parameters.  It  reinforces  the  conjecture  of  the 
resonance  theory. 
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I.  Introduction 

The  blood  pressure  is  generally  thought  to  originate 
from  resistance  of  blood  flowing  through  vessels,  which  can 
be  described  conceptually  by  Poiseuille’s  equation: 
P(pressure)  =  Qlflow^Rjresistance).  Circulation  theories 
deduced  from  the  Navier-Stokes  equation,  what  we  call 
"flow  theories”,  is  a  combination  of  two  forces:  the  inertial 
term  and  the  viscous  term  [1,  2].  In  electric  analogue,  a 
segment  of  vessel  can  be  simulated  by  series  of  an 
inductance  (L,  inertial  term)  and  a  resistance  (R,  viscous 
term). 

In  the  conjecture  of  the  resonance  theory,  the  arterial 
system  is  treated  as  a  pressure -wave-transmitting  system, 
the  pulse  pressure  arriving  at  the  terminal  will  drive  the 
blood  flux  into  the  microcirculation  [3-10].  Therefore  the 
radial  dilatation  of  the  vessel  wall  was  emphasized,  and  the 
pressure  wave  is  transmitted  in  the  form  of  “moving 
windkessel”  [3-5].  We  have  derived  an  analytic  pressure 
wave  transmission  equation  in  a  long  elastic  tube  [4,  5], 
From  this  equation,  we  also  developed  a  semi -empirical 
procedure  to  describe  the  pressure  distribution  in  a  complex 
simulated  model  composed  of  a  main  tube  and  attached 
organs  with  different  sizes  and  at  different  sites  [3].  We  got 
good  result  in  describing  the  pressure  distribution  in  these 
simulating  experiments. 

Based  on  these  previous  works,  now  we  try  to  apply 
these  fitting  methods  to  hemodynamics  in  vivo.  Since  the 
kidneys  occupies  the  largest  cardiac  output  in  all  organs,  we 
treat  the  arterial  system  as  a  system  composed  of  main 


artery  and  two  kidneys,  as  shown  in  Fig.  1  (details  see 
THEORY  section).  When  we  perform  renal  ligation,  one 
kidney  won'  t  participate  in  the  coupled  resonance  with  the 
main  artery,  thus  lead  to  the  alteration  of  blood  pressure 
waveform.  Similarly,  when  we  bent  the  aorta  transversely, 
the  elastic  property  of  the  main  artery  will  be  changed, 
therefore  the  coupled  resonance  between  the  main  artery 
and  organs  will  also  be  altered,  also  lead  to  changes  of 
blood  pressure  waveform. 

Here  we  performed  these  two  experiments  to  the  rats. 
We  observed  the  alteration  of  blood  pressure  waveform  and 
applied  semi-empirical  fitting  procedure  deduced  from  the 
resonance  theory  to  the  experimental  results.  It  can  serve  as 


Figure  1:  arterial  system  composed  of  main  artery  and  two  kidneys 
II.  THEORY 

Part  1 :  a  segment  of  aorta  in  the  resonance  theory 

In  the  conjecture  of  the  resonance  theory,  the  blood 
pressure  is  generated  by  the  momentum  of  the  injecting 
blood  from  the  heartbeat  colliding  on  the  vessel  wall,  it 
travels  along  the  aorta  in  the  form  of  transverse  wave.  The 
wave  propagation  is  like  a  “moving  windkessel”,  an  elastic 
chamber  moving  down  the  artery  with  finite  speed.  We  can 
derive  the  wave  propagation  equation  based  on  the  radial 
movement  for  a  small  piece  of  wall  [4,  5]: 


while  £  =£  +£e  ,  \Ejl=k\ 


E  :  radial  Young's  modulus  of  the  wall  material 

£e :  circumferential  Young's  modulus  of  the  wall  material 

\i  =  pwh  +  pfhf,  p„,  :  density  of  the  tube  material,  pf: 

density  of  the  adherent  fluid  which  moves  together  with  the 
wall  of  thickness,  tif ,  h  :  tube  thickness 
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:  the  inner  radius 


c  _  d{nr„~ )  :  the  arterial  compliance 

4  dP 

E  :  the  shear  modulus  of  the  tube  material 
g,  :  the  viscosity  constant  for  radial  movement  of  the  wall 
and  the  adherent  fluid 
We  will  have: 

h  4^  +  “T  +  »P  =  ~V-V~k2p 
dr  dr 

ik  =  — 1,L  =  wavelength)  (E2) 

X 

at  some  particular  position  (z  =  Z  )• 

The  term  j xvjk  2P  is  the  external  force  p  .  This  term 
includes  the  heart  pumping  and  all  the  other  effects  from  the 
vascular  bed. 

Let 

L=  p. ,  equivalent  inertance  from  a  segment  of  aorta 
/?=£,  ,  equivalent  resistance  from  a  segment  of  aorta 
1/C=  pV  j  .equivalent  compliance  of  a  segment  of  aorta 
It  becomes: 

l£L+r1L+L=L  (E3) 

dr2  dt  C  C 

In  electric  analogue,  the  left  side  of  (E3)  can  be 
represented  by  series  of  an  L  (inductance),  an  R  (resistance), 
and  a  C  (capacitance). 


Part  2:  composite  system  of  aorta  and  two  symmetric 
kidneys  in  resonance  theory 

To  take  the  kidneys  into  the  system,  we  propose  a 
pressure  coupling  mechanism  as  shown  in  Fig.  1.  The  renal 
arteries  can  be  considered  as  a  channel  of  pressure  balance, 
and  this  pressure  difference  between  the  aorta  and  the 
kidneys  is  acting  through  the  pressure  coupling  constants 


We  may  thus  construct  the  coupled  oscillation  equations 
for  the  kidney  system,  which  comprises  a  segment  of  aorta, 
two  assumed  symmetric  kidneys,  and  two  renal  arteries 
serve  to  connect  these  three  elements: 

-(Pt-P2)=-F 


,  d~P.  „  dP.  e  e 

L  — -L  +  ft,  — L  +  _L  P.  +  2- 
dr  dt  C,  C, 


C, 


,  d2P \  „  dP.  e,  „  e.,  „ 

L, — J-+R,—L+-LP^  +  -U(P.-PI)  =  0 
-  dt2  2  dt  C,  2  C,  2  1 


(E4) 

(E5) 


E  :  the  applied  force  from  the  heart  that  is  transmitted  to  this 


segment  of  aorta 

p  ,  pe-  the  pressure  in  the  segment  of  aorta  and  in  the  vascular  bed 


of  the  kidneys 

L  ,  R]  ,c  :  equivalent  inertance,  resistance,  and  amipliance  from 
the  segment  of  aorta 

L2 .  R, ,  c2  equivalent  inertance,  resistance,  and  compliance  from 
the  vascular  bed  of  the  kidneys 

e  :  the  equivalent  elastic  coefficient  of  the  whole  vascular  wall  as 


one  unit 

e.  (P.  -P.)  '■  the  coupling  force  between  the  kidneys  and  the 


segment  of  aorta,  p.  is  their  coupling  constant 


Figure  2:  Electric  analog  for  the  circulatory  system:  the  main  artery,  one 
organ,  and  the  vein. 

Similarly,  in  electric  analogue  as  shown  in  Fig. 2,  the  left 
side  of  (E4,  the  aorta)  can  be  represented  by  a  series  of  an 
RLCj,  while  the  left  side  of  (E5,  the  kidney)  represented  by 
another  series  of  RLC2  parallel  to  the  aorta.  The  return  line 
is  the  vein. 


Part  3:  Solving  the  resonance  equations  for  renal  ligation 
experiment 

Vi0!=A.  •  -1.2 

L  LC 


set  l  =  r- 


a  ~  = 


2  _  e12 


L\C\ 


’  a,"  = 


2  _  g12 


If  the  constants  in  equation  (E4)  and  (E5)  are  not 
pressure -dependent,  the  equations  become  linear.  We  may 

S6t  Fex,  =  '£FUn(0)ej"'"  P,  =  £ P, ( jnm )e i"m 

i  =  l,2,  j=4^l 

Let  m  =ru 0 


The  solutions  of  the  above  equations  are: 
F(- co„2  +  ik2  ton  +v  202  +a22)  (E7.1) 


P, ,  =- 


Fa , 


P.  2  =Pn  3=- 


with  F  _  Fex 


(E8) 


LC, 


If  one  kidney  is  completely  isolated  from  this  system, 
the  solution  will  become: 

p  '  F( — co„_  +  /LuL  +v20  +012  ) 

( — co„-  +  zk|(0n  "1-Vjq  -h cx 22 — )( — u)„- ?L co„  -cVjq-  +cx,  ) — a„-a2 

(E7.2) 

We  may  take  the  ratio  of  pnl  and  p  y  to  eliminate  the 

F  term  and  to  predict  the  effects  of  renal  ligation  on  the 
pressure  pulses. 


Part  4:  Solving  the  resonance  equations  for  aorta  bending 
experiment 

When  we  bent  the  aorta,  the  properties  of  the  kidneys 
was  assumed  unchanged,  while  those  of  the  aorta  ( ^  ,v10,  a12) 

were  shifted  to  (^',  v ,  a12')-  The  solution  for  the  blood 
pressure  wave  before  aorta  bending  was  (E7.1),  while  that 
after  can  be  represented  by: 


p  II  f(— <On  +  +V20  +  O',  ) 

"  [-to,,2  +!&j'co„  +  (v]0')2  +(a12')2](-a)„2  +0^(0,  +v202  +a22)-fa12l)2a22 

(E7.3) 

Similarly,  we  may  take  the  ratio  of  p  and  p  "  to 

eliminate  the  F  term  and  to  predict  the  effects  of  aorta 
bending  on  the  pressure  pulses. 

m.  Materials  and  methods 

1 .  Animal  preparation 

7  male  Wistar  rats  weighted  265.71±50.88  gw  (mean  ± 
SD)  were  anesthetized  with  urethane  (1.2mg  g"1  body 
weight).  Abdominal  aortic  pressure  was  measured  through  a 
PE  (polyethylene)  tube  (2F,  outer  diameter  0.6mm)  filled 
with  liquids  (0.9  %  NaCl  and  0.3  %  heparin)  inserted  from 
the  caudate  artery.  The  head  of  the  tube  was  then  moved  to 
the  level  of  the  abdominal  aorta.  Pressure  pulses  were 
recorded  by  a  catheter-tip  pressure  transducer  (P10EZ 
transducer,  RP-1500  Narco-Bio-System).  The  frequency 
response  of  this  assembly  system  was  uniform  up  to  60  Hz. 
The  signal  was  connected  to  a  simultaneous 
sample  -and-hold  panel  AX753  and  then  to  a  16-bit  AD 
converter  AX5621.  These  cards  converted  signals  to  digital 
data  with  sampling  frequency  1500  Hz,  and  sent  them  to  a 
personal  computer  for  further  analysis  [8-10]. 

In  aorta  bending  experiment,  we  opened  the  ventral  side 
of  the  rat,  separated  the  inferior  mesenteric  artery  and  tied 
prepared  silk  suture  to  it.  The  other  end  of  silk  suture  was 
tied  to  a  stereotaxic  to  provide  a  precise  transverse 
displacement.  When  we  bent  the  aorta,  the  renal  artery  was 
not  moved  at  all.  Blood  pressure  signal  was  recorded  before 
and  immediately  after  the  bending.  Right  after  the  bending, 
the  blood  pressure  was  recorded  again  to  check  if  it  was  the 
same  as  that  before. 

In  the  renal  ligation  experiment,  we  separated  the  left 
renal  artery  off  the  fat  and  then  prepared  the  silk  sutures  to 
go  around  the  artery.  We  then  threaded  both  ends  of  the  silk 
suture  through  a  segment  of  PE  tube.  When  we  performed 
the  ligation,  we  lifted  the  silk  sutures  and  confined  the  PE 
tube  at  the  same  position  to  keep  the  renal  artery  from 
moving  and  thus  restricted  the  aorta  from  bending.  The 
transverse  displacement  of  the  renal  artery  was  far  less  than 
1  mm.  Recording  procedure  was  the  same  as  that  in  aorta 
bending  experiment. 

After  the  surgical  operation,  we  waited  30  minutes  for 
the  blood  pressure  wave  to  return  to  normal  to  start  our 
experiments.  At  least  15  minutes  were  waited  between  each 
step  to  avoid  interference  from  the  last  experiment.  Between 
the  experiments,  the  exposed  tissues  were  kept  moist  with 
physiological  saline  and  the  environment  was  kept  at  32°C 
by  a  heat  radiator. 

2.  Data  analysis 

In  each  recording,  we  have  10  pulses.  Each  pulse  was 
picked  by  two  relatively  lowest  points  and  transformed  into 
frequency  spectra  by  Fourier  transform.  To  eliminate  the 
effects  of  time  -variant  F,  amplitude  of  each  harmonic  was 


normalized  by  its  mean  level  (DC  component)  and  named 
harmonic  proportion.  Because  the  amplitude  decreased 
rapidly  with  the  harmonic  numbers,  we  focused  only  on  the 
first  six  harmonics  [8-10]. 


Figure  3:  Harmonic  variations  of  pressure  wave  by  renal  ligating.  Upper: 
amplitude.  Lower:  phase  angle.  Solid  line:  measured  data.  Bars  are  SD. 
Dash  line:  fitting  by  the  resonance  theory. 

The  fitting  procedure  was  similar  with  our  previous 
work  [3].  We  first  compare  the  measured  and  theoretic 
spectrum  in  renal  ligation  experiment  and  apply  the 
least-mean-square -error  criteria  to  both  amplitude  and  phase 
angle  of  first  4  harmonics  to  get  6  parameters:^,  £2(=l), 
v10>  v20(=v30),  a12(=  a13)  and  «,(=«,)■  Similarly,  we  then 
turn  to  the  result  of  aorta  bending  experiment  and  get 
another  9  parameters:  k.t '  (5mm) »  k^' { 6mm) >  v10'(4mm)> 

V 10 ' (5mm),  V  l0'  (6mm),  cc12'(4mm),  ot12'  (5mm)  and  ci12'(6mm)  (inside 

the  parenthesis  is  the  transverse  bending  of  the  aorta). 

In  all  these  equations,  the  F  term  from  the  heart  is  a 
problem,  because  it  is  an  individual-dependent  variable. 
Different  rat  will  have  different  heartbeat,  even  for  the  same 
rat,  this  term  can  be  time -dependent.  However,  Tsecond 
bending  is  quick  enough  to  avoid  any  reflex  reaction  and  F 
can  be  treated  as  a  constant.  This  is  confirmed  by  carefully 
checking  the  pressure  pulses  before  and  after  the  bending  to 
be  the  same.  In  the  theoretical  treatment,  under  the  premise 
that  the  intervention  does  not  change  F  ,  we  may  take  the 
ratio  of  pm  and  pm  (or  pm  )  to  eliminate  the  F  term, 
thus  to  predict  the  effects  of  renal  ligation  and  aorta  bending 
on  the  pressure  pulses. 

When  we  compare  the  first  6  harmonics,  both  with 
amplitude  and  phase  angle,  we  have  12  equations. 
Considering  4  conditions  of  renal  ligation  and  bending  the 
aorta  for  4,  5,  6  mm,  we  have  totally  12x4=48  equations, 
while  there  are  only  15  independent  fitting  parameters. 

IV.  RESULTS 

Comparison  between  experimental  data  and  curve  fitting 
by  the  resonance  theory  in  renal  ligation  experiment  is 
shown  in  Fig.  3,  and  that  for  aorta  bending  is  shown  in  Fig. 


4.  The  fitting  parameters  in  Fig.  3  are  fo  =  6.0 Hz  ,  <b0  =  2n/0 , 
^  =70,  fc,  =k3  =  60>vlo  =  co0 ’V20  =v30=  2.6co0,  a122  =a132  =  3co02 , 
and  a22=a32  =  1.5co02  ( to0 ,  v,  k  and  a  are  in  radians). 

The  fitting  parameters  for  Fig.  4  are  the  same  as  those  in  Fig. 
3  except: 

(A)  v10'=1.6vlo>  al2'=  0.83al2>  A:, '=  1.03fe, i 

(B)  v10'=Zlv10>  a12’=0.6a12,  ^'=1.05^; 

(C)  v10’=2.27v10,  a12'=  0.45a12 ,  ^'=1.1*,- 

V.  DISCUSSION  &  CONCLUSION 

In  the  fitting  results  (Fig.4),  the  three  changing  parameters: 
v  |o' ,  (x  ’  at  transverse  displacement  of  4,  5,  6  mm  show  rational 
trends,  /t, '  increased  from  1.03 /t,,  1 .05  k{ ,  to  1.1^,  v10'  from 
1.6Vl0,2.1Vl0,to2.27v10,and  «12'  changed  from  0.83  0t|2 , 0.6a,2 
to  0.45  K| , .  It  is  consistent  with  intuitive  thought  that  when  the 

aorta  was  bent  more,  these  parameters  deviated  more  from  their 
unbent  value  accordingly.  Besides,  our  fitting  parameters  are  not 
chosen  arbitrarily  but  have  their  physical  meaning,  therefore  we 
call  them  '‘semi-empirical”  ones  [3-5], 

Moreover,  48  linear  equations  can  generally  determine  48 
parameters.  However  in  our  results,  we  can  use  only  15 
independent  parameters  to  fit  totally  48  equations.  We  can  see  that 
the  48  equations  are  so  highly  cbpendent  that  these  48  parameters 
can  be  reduced  into  15  independent  ones  as  suggested  in  our 
description.  It  illustrates  that  even  if  there  are  so  many  complex 
structures  in  the  arterial  system,  we  can  still  get  good  description 
for  the  blood  pressure  distribution  by  dividing  the  arterial  system 
into  sub-units  and  describing  with  few  elastic  parameters  as  in  the 
resonance  theory. 

This  study  serves  as  the  toughest  test  for  our  semi-empirical 
approach  in  vivo ,  and  the  result  elucidates  the  conjecture  of 
the  resonance  theory. 

Figure  4:  Harmonic  variations  of  pressure  wave  when  the  aorta  was  bent 
for  a  transverse  displacement  of  (A)  4  mm;  (B)  5  mm;  (C)  6  mm.  Upper: 
amplitude.  Lower:  phase  angle.  Solid  line:  measured  data.  Bars  are  SD. 
Dash  line:  fitting  by  the  resonance  theory. 
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